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SUMMARY 


A post-processor system for the data reduction and post analysis of NASTRAN results is 
described. NASTRAN analysis results are scanned to determine maximum and minimum displace- 
ments, forces and stresses. Allowables and margins of safety are computed, and in the case of 
multiple loading conditions, envelopes for displacements, forces, stresses and margins of safety are 
also produced for specified element sets. Graphical plots of the reduced or the regular NASTRAN 
results may be obtained superimposed either on a developed fuselage strip or on a projection of 
any specified part of the finite element model. The use of the data reduction, post analysis and 
graphical plotting capabilities provide the analyst with a fast and convenient tool for the study of 
NASTRAN analysis results and their presentation for project documentation. 


INTRODUCTION 


The high computational speed and large storage capacity of modern computers have enabled 
the analysis of large and complex structures subjected to multiple loadings and various boundary 
conditions. As a consequence, the structural engineer is bound to devote a substantial portion of 
his time in scanning and interpreting a large amount of output data, a process which is time con- 
suming, error prone and hence inefficient. The way to alleviate this problem is to automate, 
wherever possible, the scanning and interpretation of the results, and to give the analyst the option 
to reduce the amount of computer output according to his engineering requirements. For further 
efficiency, this should be done in conjunction with graphical display of the reduced finite element 
analysis results. Additional help for the analyst can be obtained through the automation of certain 
standard post analysis procedures such as the computations of allowable stresses and margins of 
safety for the structural elements used in the finite element analysis. 

The need for these capabilities has been felt for some time and as a result a large number 
of post-processors have been developed. However, most of these satellite programs are generally 
limited in scope, usually serving one particular purpose only, such as plotting (reference 1), or 
scanning (reference 2). At Israel Aircraft Industries an attempt has been made to integrate data 
reduction, post analysis, and graphic visualization into one package which has been especially 
devised to be engineering oriented. 

This paper describes two interrelated software modules for the post processing, post analysis 
and graphical presentation of NASTRAN analysis results. Emphasis is put on the fact, that both 


modules are easy to apply, demanding a minimum of information from the user, yet leaving him 
with enough options to answer his design oriented requirements. The first of the two modules, 
NASDAT (reference 3) performs the post-processing and post analysis of the NASTRAN results; 
the second module, GRAS, plots any of the regular or reduced results in a graphical form. 


DESCRIPTION OF MAIN MODULES 
Description Of The NASDAT Module 

The execution of the NASDAT module, illustrated in figure 1, may be subdivided into three 
main phases. 

a) In phase one, regular NASTRAN output data blocks are sorted and stored in a compact 
form on a disk file for subsequent use. This file may be saved after termination of 
NASDAT execution. Thus, a restart file, which contains all necessary I/O information 
of the NASTRAN analysis, is created to be used in phases 2 and 3. 

b) In phase two, the packed and sorted NASTRAN output is scanned, and efficient 
data reduction performed in accordance with the user’s requirements. 

c) Finally, in phase three, post analysis of results including computations of allowables 
and margins of safety are performed. 

NASDAT may be run either in conjunction with NASTRAN as a single job, or as a separate run 
following a NASTRAN analysis for which the necessary output data blocks have been saved. 
NASDAT may be restarted as many times as required. 

NASDAT Capabilities 

Presently, using NASDAT capabilities, the user may obtain: 

* Maxima and minima values 

* Envelopes of displacements, forces, stresses and/or margins of safety 

* Output in preferred sequence and format 

* Computations of allowables and margins of safety 

The following is a brief description of the different options mentioned above: 
a) Maximum and Minimum Values 

The user may define a given set of elements and/or gridpoints from which the program 
will single out and identify those elements and grid points bearing extreme values of 
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force, stress and displacement; the magnitude of these extrema is also produced. 

In addition, it is possible to specify the upper and lower bounds beyond which the 
search for maximum and minimum values is to begin. 

b) Envelopes of Displacements, Forces, Stresses and Margins of Safety 

When analyzing several loading cases or checking various boundary condition con- 
figurations, the user has the option of obtaining envelopes of displacements, forces, 
stresses and margins of safety for a specified set of grid points and elements. 

c) Output In A Preferred Sequence And Format 

In some cases the user may require the results in a certain region of interest to be 
printed out separately in a specified sequence. In addition, if upper and lower bounds 
on the output values are also defined, then only those elements of the set with results 
above these bounds are printed. 

d) Computations Of Allowables And Margins Of Safety 

The margins of safety of any Rod, Shear Panel or Membrane element can be computed 
using allowables either defined by the user or automatically computed by the program. 
For the latter option the user may either define the section properties via manual input 
or allow the program to retrieve the geometrical and mechanical properties of these 
sections from tables generated in the NASTRAN analysis. 

The NASD AT Input 

In order to generate the main input for NASDAT the following DMAP statements for 
static analysis, must be included into the Executive Control Deck of the NASTRAN deck: 

1. ALTER 121 

2. OUTPUT2 

3. OUTPUT2 

4. ENDALTER 

Statement 2 saves the necessary geometrical and mechanical properties of the structure on 
NASTRAN Fortran file UTl and Statement 3 saves all the analysis results. 

As can be noted, the main input to NASDAT comes from the above mentioned data blocks 
stored on file. However, additional input is required to specify the user's request for data reduction, 
post analysis and graphical presentation of results. This is prepared in a form similar to the 
NASTRAN input (reference 4). The NASDAT data deck consists therefore of an Executive Con- 
trol Deck, a Case Control Deck and a Bulk Data Deck. 


CASECC, EST, MPT, GPL,/C,N,- 1 /C,N, 1 1 /C,N,UTAPE $ 
OUGV 1 , OQG 1 , OEF 1 , OES 1 ,/C,N,0/C,N, 1 1 $ 


- r- -» 

:>oi 


The Executive Control Deck contains that information required by NASDAT to decide if 
the NASTRAN output files are to be sorted and/or checkpointed for restart purposes. 

The Case Control Deck consists only of SET cards needed to define the element and grid 
sets referred to in the Bulk Data Deck, and TITLE cards which produce output labels for the 
preferred sequence option. 

The Bulk Data Deck contains the main information required for the data reduction and 
post analysis. The format and order of input is similar to that of NASTRAN, i.e. data is input 
in 8 column fields with freedom to arrange the input cards in any sequence. Facsimiles of input 
cards for the data reduction and definition of lower and upper bounds are shown for Bar and 
Shear elements in figures 2, 3 and 4 respectively. Facsimiles of input cards for the automatic 
computation of margins of safety for Rod and Shear elements are shown in figures 5 and 6 
respectively. 


Description Of The GRAS Module 

This is the off-line graphical postprocessing module, which is used to display graphically 
either the reduced results produced by NASDAT or the regular NASTRAN output. Input for 
the GRAS module is generated by NASDAT. 

At present, GRAS is run in a batch mode environment and uses a Calcomp drum plotter. 
The plotting capabilities can take one of the following forms: 

a) Developed Fuselage Strip 

Using this option, a segment of the fuselage of an aircraft is developed into a flat 
strip upon which the values of the reduced or regular results are plotted. 

b) Bulkhead Plots 

Using this option graphical displays of bending moments, bending stresses and trans- 
verse shear forces in Bar elements may be obtained. 

c) Projection Plots 

The user may obtain a graphical display of the projection of any specified part of 
the structure on any one of the three principal planes. The stress values are then 
plotted on this projection. Figure 7 shows stresses in Rods plotted on a projected 
view of an aircraft delta wing in the XY plane. 
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NUMERICAL APPLICATIONS 


The NASDAT and GRAS modules were first utilized in the data reduction and post 
analysis of the results obtained from the NASTRAN analysis of the lAI Westwind 1124 aft 
fuselage (reference 5). The structure was idealized using Shear Panels to represent the skins, Rod 
elements to represent the stringers and Bar elements offset from the skins to represent the bulk- 
heads. The complete mathematical model, shown in figure 8, consisted of 460 BAR elements, 

1670 Ro 6 e/ements, 1030 SHEAR elements and 1250 GRID points representing a total of 3650 
unconstrained degrees of freedom. The structure was subjected to 18 loading conditions. Because 
of the size of the structure a substructuring technique (reference 6) was employed, subdividing 
the aft-fuselage into 5 substructures. 

Based on engineering considerations related to the structural configuration and loading 
conditions the analyst has asked for (presented here as an example) information relative to extreme 
values of average shear stresses, data relative to critical regions affected by shear buckling of skin 
panels, and margins of safety of Rod elements. In order to get the full picture of the internal load 
distribution in a bulkhead of interest, the preferred sequence option has been used to output re- 
sults in this region. In addition, graphical capabilities have also been employed to display moments, 
shear stresses and axial stresses. Finally, a plot to represent the envelope of shear stresses has been 
required. 

A listing of the complete input required to perform this type of data reduction and post 
analysis for substructure 1 (figure 9), is given in figure 10. Based on the user’s request, the fol- 
lowing data was obtained: 

(i) Twenty extreme absolute values for average stresses in Shear Panels and axial stresses 
in Rod elements (figure 1 1 ). 

(ii) Buckling stresses, diagonal tension factors and related margins of safety for Shear 
Panels located in the left hand side of the fuselage structure. 

Since the actual skin panels between bulkheads were represented in the idealization 
as a mesh of NASTRAN Shear Panels, it should be noted that the aspect ratio 
(a/b) referred to the actual skin panels and not to the dimensions of an idealized 
element. Therefore, in this case, the user had to input the aspect ratio manually. The 
critical stresses, diagonal tension factors and margins of safety, shown in figure 12 
were computed using criteria defined in references 7 and 8. 

(iii) Euler buckling stresses, ultimate tensile stresses and related margins of safety for Rod 
elements located in the left hand side of the fuselage structure. 

In this case the necessary geometric material and cross section properties of the rod 
elements were retrieved from the NASTRAN tables using the automatic procedure 
(via the FLAG=2 option). Computer' output is shown in figure 13. 

Note that this output is presented for illustration only. Usually, the user will request 




only a print-out for those elements having negative margins of safety. 

(iv) A table of the internal load distribution for a bulkhead as requested via the preferred 
sequence option. Computer output is shown in figure 14. 

(v) Graphical display needed for project documentation including (a) plot of bending 
stresses for a bulkhead as shown in figure 15; (b) presentation of envelopes for 
maximum and minimum shear stresses as shown in figure 16 and 17 respectively, 

(c) presentation of stresses in Rod elements and Shear Panels as shown in figures 

1 8 and 1 9 respectively. 

Since the original height of the Calcomp plots was 29 inches (73.6 cm) only a part of 
the plot (the left hand side) is presented. 

The total amount of computer output for the complete NASTRAN finite element analysis 
consisted of about 5000 pages. Scanning, hand reduction, post analysis and drawing of results 
would have required several man months of engineering work. In contrast, meaningful engineering- 
oriented output and data reduction as presented above could be accomplished in the span of a few 
days. 


CURRENT DEVELOPMENTS 

Current developments of NASDAT include the expansion of the options which deal with 
computations of the allowables. These will incorporate bending elements as well as stiffened panels. 

The NASDAT/GRAS modules are presently converted to operate in a time-sharing mode 
using low cost graphic display terminals, and together with NASTRAN will be included as modules 
in Israel Aircraft Industries’ ISSAS, (Interfaced Software System for Sizing and Analysis of air- 
craft Structures) (reference 9). 


CONCLUSION 

Scanning of analysis results and computations of allowables and margins of safety, which 
are especially error prone and time consuming, may be performed with the aid of NASDAT and 
results grap.hically displayed by GRAS in an automated fashion, thus freeing the analyst to devote 
a larger portion of his time to engineering decision-making based upon results which are obtained 
in an organized form. 

The use of the presented data reduction, post analysis and graphical plotting capabilities 
provide the analyst with an efficient and convenient tool for the study of NASTRAN analysis 
results and their presentation for project documentation. 
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FIGURE 1 : FLOW DIAGRAM FOR THE NASDAT MODULE 
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Bulk Data Card Data reduction for CBAR elements l/MINMAX Nuirtjer of minimum or maximum extreme stress or 

Description: Defines the option and the element set for force values required (Integer > 0 or blank) 

which data reduction is requested. Generate input for the plot modules 

"YES" - generate 

Format and Example 

"NO" - do not generate (Default) 
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FIGURE 2: NASDAT BULK DATA CARD BAR 
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FIGURE 3: NASDAT BULK DATA CARD BARMIN 
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FIGURE 4: NASDAT BULK DATA CARD SHEAR 






Description Defines the properties of CROD and CONROD 
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FIGURE 5: NASDAT BULK DATA CARD ALROD 




Bulk Data Card AISHEAR Computations of Allowables and related SUBCASE Single Subcase identification number for which 

margins of safety for Shear elements. data reduction is requested (Integer > 0) 
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FIGURE 6: NASDAT BULK DATA CARD ALSHEAR 
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FIGURE 8: FINITE ELEMENT MODEL OF WESTWIND 1124 

AFT FUSELAGE (COMPLETE MODEL) 




(CSTWIND-SUKTRUCTURT 1 


UNDr> QRMEO SHW=E 

FIGURE 9: FINITE ELEMENT MODEL OF WESTWIND 1124 

AFT FUSELAGE (SUBSTRUCTURE 1) 
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NASOAT EXECUTIVE CONTROL DECK ECHO 


10 POST* PROCESSOR 
SORT YES 
CHKPNT YES 
CEND 


.MASQAT CASE CONTROL DECK ECHO 


SET 1*1121 thru 1162 Except UO*U<»6#ll5l*ll52.1160tU70*n74tllT5t 
1179 thru 1186*1199 THRU 1206*1695*1557 
SET 2*1 163* I 166* 1151 *1152* 1160*1168* 1169* 1172* 1173, I 176* 1177*1178* 

1287, 1501 *1507* 1521 *1527*1561 *1565* 1569 
SET 3«116S*1166*1167*ll&5 THRU 1190 
SET 6*1191 thru 1198*1553 

sIt Hmis’tmru 3118*3181*3182*3189*3190.3199 THRU 3212 EXCEPT 3205*3206* 
3557*3561 

SET 0*3lll* 3122*3129*3130*3166*3167*3171 *3172*3176*3177*3187*3188*3195* 
3196*3569*3553.3557,3561 

SET 9«3l23 Thru 3128*3135*3136*3166*3168*3169*3691*3565 
SET 10*3138*3166*3191 THRU 3196 

SET 11*3139 THRU 3151 EXCEPT 3165*3166*3503*3517*3523*3529*3539 

SET 12*3156*3162 

SET 13*3197*3198*3205*3206 

SET 16*3152 THRU 3162*3173 Thru 3175*3178 THRU 3180 
SET 15*3153 thru 3186 

titlEhTrahe at station 360 - substructure 1 - 
SET 25*2 THRU 11 
SET 21*57 T^RU 86 
8E61N BULK 


N A S 0 A T 


SORTED BULK 


DATA ECHO 


CARO 




COUNT 




1 

ALROO 

1 

0*132 

2 

ALROO 

2 

0*167 

3 

ALROO 

3 

0.68B 

6 

ALROO 

6 

0*652 

5 

ALSHEAR 

S 

6*0 

6 

ALSHEAR 

6 

3*8 

7 

ALSHEAR 

7 

3*6 

8 

ALSHEAR 

8 

3*0 

9 

ALSHEAR 

9 

2*0 

10 

ALSHEAR 

10 

1*5 

11 

ALSHEAR 

11 

1*0 

12 

ALSHEAR 

12 

1*0 

13 

ALSHEAR 

13 

5*0 

16 

ALSHEAR 

16 

2*5 

IS 

ALSHEAR 

15 

1*5 

16 

8AR 

21 

15678 

17 

ROD 

ALL 

2 

18 

ROO 

ALL 

2 

19 

ROOM IN 
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20 

SHEAR 

all 


21 

SHEAR 

ALL 


22 

SHEAR 

ACL 


23 

enmata 




2 

2 

2 

2 


0*06 

10*5E«6 

0*3 


1 

0*06 

10*5E*6 

0*3 
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0*06 

10*SE«6 

0*3 


1 

0*06 
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0*3 
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0*06 
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0*06 

10.5E«6 
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0*06 

10*5E«6 
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1 

0*06 

10.5E^6 

0*3 


1 

0*072 

10*5E*6 

0*3 


1 

0*872 

10.5E*6 

0*3 
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FIGURE 10: TYPICAL INPUT DATA DECK FOR NASDAT 
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THE f0tL0WlN6 tI$T P«EsCNT$ T»<C EXTRCNC STtCSSCS IN € S H E A R CLEMENTS FOR S U • C A S E 
20 EXTREME VALUES ARE REOUESTEO - 10 VALUES ARE FOUND TO 8E POSITIVEt 10 


ARE negative 


MAXIMUM IFOSITIVEI STRESSES 


CLEMENT 

AVG 

10, 

STRESS 

3151 

1,1645E«04 

3196 

1«0367E^04 

3149 

9«9385E«03 

3194 

9.5829E^03 

3166 

8.9071E<»03 


element avg 

ID. STRESS 

3150 1«0372E^04 

3156 1*0170E«04 

3195 9«8036E*03 

3193 6.9942E^03 

3168 8«5570E«03 


MINIMUM (NEGATIVE) STRESSES 


CLEMENT AVG 
ID. stress 
3049 •9.4084E«03 
3073 -7*7787E^03 
3040 •7«6418C^03 
3056 -7,3720E^03 
3036 -7.3221E^03 


CLEMENT AVO 

SO, STRESS 

3066 •6«4594E^03 
3048 •7,6786E«03 
3060 -7.4885C«03 
3032 -7.3580E«03 
3047 -7,3046E«03 


the following list PRESENTS TME EXTREME AXIAL STRESSES IN C R 0 0 ELEMENTS FOR SU8CASE 4 
20 extreme values are REQUEST - 20 VALUES ARE FOUND TO 8E POSITIVEt -0 ARE NEGATIVE 






maximum (POSITIVE) 

stresses 






element 

axial 

CLCHENT 

AXIAL 

ID, 

STRESS 

ID, 

STRESS 

1139 

4,9660C^04 

1097 

4,8070C«04 

1138 

4,7072E^04 

1089 

4,6432E^04 

1113 

«.63S3E*04 

1090 

4,6018E^04 

1121 

4,5952E«04 

1129 

4,5700C^04 

1091 

4,5455E^04 

1114 

4,5029E«04 

1122 

4,4157C^04 

1134 

4,4011E^04 

1115 

4,3012E«04 

1123 

4,2808C«04 

1105 

4,2489C«04 

1133 

4,1843E«04 

1102 

4,17S7E^04 

1140 

4,1635E404 

1124 

4,0831C«04 

1U5 

4,0641C^04 


FIGURE 11: MAXIMUM AND MINIMUM STRESSES IN 
SHEAR PANELS AND ROD ELEMENTS 
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FIGURE 12: AVERAGE STRESSES, BUCKLING STRESSES, 

DIAGONAL TENSION FACTORS AND MARGINS 
OF SAFETY FOR SHEAR PANELS 
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FIGURE 13: AXIAL STRESSES, EULER/ULTIMATE STRESSES 
AND MARGINS OF SAFETY FOR ROD ELEMENTS 
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FIGURE 14: INTERNAL LOAD DISTRIBUTION IN A BULKHEAD 


STREbSEb IN BAR ELEMENTS 
BULKHEAD AT STATION ^ 340.00 


FIGURE 15: BENDING STRESSES DIAGRAM FOR A BULKHEAD 



NESTWIND SUBSTRUCTURE 1 
ENVELOPE FOR MAXIMUM 
SHEAR STRESS IN PANELS 


FIGURE 16: POSITIVE SHEAR STRESSES ENVELOPE DISPLAYED 

ON A DEVELOPED STRIP OF THE FUSELAGE 


WESTWIND SUBSTRUCTURE 1 
ENVELOPE FOR MINIMUM 
SHEAR STRESS IN PANELS 


FIGURE 17: 


NEGATIVE shear STRESSES ENVELOPE DISPLAYED 
ON A DEVELOPED STRIP OF THE FUSELAGE 




WESTWIND SUBSTRUCTURE h 

SUBCASE 4 -UPGUST CONDITION CONFIG. 5- 

STRESSES IN ROD ELEMENTS 

FIGURE 18: ROD STRESSES DISPLAYED ON A DEVELOPED 

STRIP OF THE FUSELAGE 


579 



WESTWIND SUBSTRUCTURE 1 

SUBCRSE 4 -UPGUST CONDITION CONFIG. 5- 

SHEflR STRESS IN PANELS 

FIGURE 19: SHEAR STRESSES DISPLAYED ON A DEVELOPED 
STRIP OF THE FUSELAGE 
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